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We consider the effect of a canonically normalised scalar field degree of freedom on the dynamics 
of gravity from small to large scales. We show that the effects of modified gravity can be completely 
captured by the time variations of the scalar field mass and its coupling to matter. This leads to 
a parameterisation of modified gravity where local constraints are easy to analyse and large scale 
scale structure effects apparent. 



PACS numbers: 

Modified gravity is an alternative to dark energy [1]. 
Plausible scenarios reproducing the accelerated cosmic 
expansion have been proposed in the past decade, and 
general parameterisations are currently sought for [2]. A 
simple possibility uses two functions v(k,a) and j(k,a) 
through the Poisson cquation[3, 4] 

-fc 2 $ = 4tt(1 + v)G N a 2 8p m 

and 

* = (1 +7)$. 

Here Gat is Newton's constant, and $ are the poten- 
tials in the Newtonian gauge: 

ds 2 = -a 2 (l + 2*)dr 2 + a 2 (l - 2<S>)dx 2 , 

where r, x are the conformal time and comoving coordi- 
nates. Because v, 7 are free functions of both time (a) and 
space (fc), such a parameterisation is very general and 
convenient for phenomcnological studies of large scale 
structure. In the linear regime, it can be shown to be cov- 
ering several important classes of modified gravity theo- 
ries, such as f(R) gravity [5]. These theories differ dras- 
tically from the v — 7 parameterisation in the non-linear 
regime where a screening of the modified gravity effects 
is present. Such a shielding effect is apparent in N-body 
simulations [6, 7]. Moreover, designed for the analysis of 
the cosmological perturbations in the linear regime, the 
v — 7 parameterisation is not applicable to other gravita- 
tional regimes such as the solar system or laboratory tests 
of gravity; in particular, it fails to capture the environ- 
mental dependence of screened modified gravity models 
that is so crucial to evade local constraints [8]. Also, it 
focuses on the phenomcnological consequences of modi- 
fied gravity and makes the physics rather obscure, e.g., 
it gives no hint about the theoretical properties of the 
extra degrees of freedom mediating the modification of 
gravity. Finally, it is not clear that a simultaneous pa- 
rameterisation of the time and spatial dependence of v, 7 
is consistent with any underlying theory. All in all, the 



v — 7 parameterisation is useful for large scale structure 
in the linear regime but cannot, in general, lead to a 
well-defined theory of modified gravity, which must be 
applicable from small scales where gravity is tested in 
the laboratory to cosmological scales where non-linear 
effects are crucial to structure formation. 

In this paper we propose a new parameterisation of 
a broad class of theories which involve a fifth force me- 
diated by a new scalar degree of freedom, such as the 
chameleon [9], dilaton [10] and symmetron [11] theories, 
and f(R) gravity. The success of these theories relies on 
certain mechanisms that suppress the fifth force in lo- 
cal, high matter-density, environments. It may seem that 
a full parameterisation of these theories should include 
not only the temporal and spatial but also the environ- 
mental dependences; in fact as we will show below, what 
we actually need is the temporal dependence of the mass 
and the coupling to matter alone, which are often simple 
power-law functions: once this is given, the behaviour of 
the theory in different regimes is completely fixed. 

The action governing the dynamics of a canonically 
normalised scalar field ^ in a scalar-tensor theory is of 
the general form 

S = J ' d'x^— g ^R- l -(V<p) 2 

+ J d 4 xv^£ m (^,^), (1) 

where g is the determinant of the metric g^, R is the 
Ricci scalar and ipm are various matter fields labeled by 
i. A key ingredient of the model is the conformal coupling 
of (f> with matter particles. More precisely, the excitations 
of each matter field ipm follow the geodesies of a metric 
g^ v which is related to the Einstein-frame metric g^ by 
the conformal rescaling 

g» v = A 2 ((j))g^ u . 

The Klein Gordon equation is modified due to the cou- 
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pling of the scalar field to matter: 

□0 = -/3T + ^, 



(2) 



where T is the trace of the energy momentum tensor T^ 1, 
and the coupling of to matter is defined by 



[3 = TOp] 



din A 

d0 ' 



(3) 



This is equivalent to the usual Klein-Gordon equation 
with the effective potential 

U off (0) = V(4>) - (A(<f>) - 1)T (4) 

In the weak-field limit with 

ds 2 = -(1 + 2$ JV )dt 2 + (1 - 2& N )dx i dx i , 

the modified geodesic equation for matter particles reads 



A 2 



-d i ($ N +]nA(<f>)). 



(5) 



where <3?at is Newton's potential. This can be interpreted 
as the motion of a particle in the effective gravitational 
potential defined as 

§ = $at + In .4(0), 

the scalar field induces a modification of gravity. 

We focus on models where the effective potential has a 
minimum 0(a) which depends on the scale factor a due 
to the time variation of the matter density. Using the 
definition of the scalar mass at the minimum of V e ff , 

2 _ 9F eff (0) 



we deduce the relation 
d 2 V 



V 



m\a)- /3 2 A(0)- 
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dft_ 
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A(0)- 
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where the couplings to matter /3 can be field dependent 
and m 2 is evaluated at the minimum of V Q g. Using the 
minimum equation 



a<p mpi 



(0) 



we find that the field value at the minimum evolves ac- 
cording to 



dt m z mpi 



(7) 



This leads to the solution for the time evolution of the 
minimum 

<K«) = zr- f S^p( a ) da + <f>c (8) 



mpi J a am 2 (a) 



where p is the density of non-rclativistic matter species, 
4> c the initial value of the scalar field and we have taken 
A((f>) w 1 as the temporal variation of fermion masses 
must be very weak (see later). If the coupling (3 is ex- 
pressed in terms of the field <fi and not the scale factor a, 
this is also equivalent to 



mpi J ai . am 2 (a) 



p{a)da. 



(9) 



Similarly the minimum equation implies that the value 
of the potential at the minimum is given by 



V{a) = V - 3 



/?(«) 



P 



am 2 (a) 



da. 



(10) 



It turns out that (10) and (8) define a parametric rep- 
resentation of the potential V(4>) and the coupling /3(0). 
This implicit dependence is valid for all the values of 4>{a) 
as a varies. In practice, this allows one to reconstruct the 
full dynamics of the model for ranging from its value 
after inflation to now. This defines the bare scalar field 
potential V(0) parametrically when /3(a) and m(a) are 
given and allows one to reconstruct the full dynamics of 
the models defined by the action (1). 

As a first example, let us consider the important case of 
a non- vanishing coupling function /3(a). Defining p = 
/3(a) = /3o<?(a) and m = mo/(a), in which a subscript o 
denotes the present-day value, we find that 



mpi 



H 2 

9/3o^mO — 9 



da- 



9(a) 



(11) 



where Hq is the Hubble constant and f2° n = /?o/3J?o m Pi- 
Let us specialise to the case where 



m(a) = uiqcl r ', /3(a) = /3o 
where r > 3/2. We find that 

0(a) - <p c _ 9/3 Q tt m0 H% 2r _ 3 2r _ 3 



mpi (2r — 3)m 



2( a ~ a ini 



(12) 



(13) 



and 



V(a) = V 



(2r - 6)m§ 
from which we deduce that 



K r " -<r b ) (14) 



U(0) = v a 



27a 2 :r^ nl om 2 l H$ 
(2r — 6)mg 



X ((1 + 0^2 (^)) a ^<^- 3 ) - 1) 

(15) 



9a?7 3 /3 fi m0 ^ 2 v mpi 



which defines inverse power law chameleon models when 
3/2 < r < 3 [16] and power law models when r > 3. In 
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TABLE I: Some f(R) gravity models, the corresponding V(<f>) 
in the Einstein frame, and the evolution of the scalar field 
mass m{a). £,n,w,A, i?» are constant parameters, with n > 
0,0 < w < 1. 



and therefore 



Ref. 



Asymptotic f(R) Equivalent V((f>) m(a)/mo 



[14,20-22] -2A+(-Sf)" 
[22] _ 2 A-£ln(£ 

[23] 



C[l-exp(-/3 K4 0)]^ a -f(-+2) 
A — Cln(«40) a -3 

Coxp(,8K 4 0) 



[l-CX P (-j3K 4 c/>)l 1 



p + ga 



the latter case and when /3o = l/y/6, these models are 
equivalent to large curvature f(R) gravity. Some choices 
of such f(R) models that lead to viable cosmologies are 
summarised in Table I, and we see that the different f(R) 
models in [14, 20-23] result in a power-law evolution of 
to with r > 3. This can be easily seen by choosing 



2A 4 

f(R) = R-^- 
m pi 



/go Rq 
n R n 



71+1 



(16) 



expanding in Rq / R for R > Rq where Rq is the curvature 
now. Using the equivalence 



and 



we find that 



h 



-2/3 o 0j i /mp 1 



m P1 Rf R - f(R) 
2 fl 



— py ( y/(n+i) 
Ro mpif Ro 



and 



V(4> R ) = k + ^fM-^^r /in+1) 



(17) 



(18) 



(19) 



(20) 



which is a power law model with n = — 2. It can also 
be obtained that the mass too is related to fn according 
to 



m<) 



mo _ _ 
H ~V (n + l)\f Ro 



(21) 



This completes the identification of the large curvature 
models with the m(a) — /3(a) parameterisation. 

As another example, consider a very different, so-called 
dilaton, model in which the coupling function /3(</>) van- 
ishes for a certain value </>* of the scalar field 4>. It is 
enough to study the dynamics in the vicinity of the field 
(/)*, where 

i Ao(6 — <£,), 



da 



from which we deduce that 



In 



9A 2 fl m0 Ho 



a 4 m 2 (a) 



(22) 



This expression will be useful to analyse gravitational 
tests of dilaton models. 

The m(a) — /3(a) parameterisation allows one to study 
all the different regimes of the models. Let us first con- 
sider the background cosmology. Constraints from the 
creation of the light elements imply that the variation of 
fermion masses during Big-Bang Nucleosynthesis (BBN) 
must be small. This implies that the scalar field must 
have settled at the minimum 4>{a) of the effective po- 
tential before the electron decoupling [12]. As long as the 
mass of the scalar field m(a) is much greater than the 
Hubble rate, m{a) 3> H(a), the minimum of the effec- 
tive potential V e s is stable. Indeed if this were not the 
case then the scalar field would receive a kick due to 
the abrupt variation of the trace of the energy momen- 
tum during the transition and induce an 0(1) change in 
the fermion masses. This is not the case if the field is 
at the minimum as the steepness of the potential, there 
m 3> H, prevents large excursions of the scalar field from 
the minimum. At the transition p ~ m\ and z e ~ 5 • 10 8 , 
the matter density at the electron decoupling is of order 
10 _5 g • cm~ 3 . We must impose that the field has settled 
at the minimum of the effective potential earlier than the 
electron decoupling for instance with z- ln i ~ 10 2 z e . At 
this particular redshift, the energy density of matter is of 
the order of lOg • cm~ 3 . Hence the initial value <f> c of the 
scalar field corresponds to the scalar field value in ordi- 
nary matter such as on the Earth, the Sun or laboratory 
test masses. Given the time evolutions of the mass m(a) 
and coupling /3(a), one can reconstruct the dynamics of 
the scalar field (j)(a) for densities ranging from cosmolog- 
ical to solar system values using Eq. (8). By the same 
token, the interaction potential can be reconstructed for 
all values of 4> (or p) of interest using Eq. (10), from the 
solar system and Earth to cosmological background: a 
tomography of modified gravity. 

The to (a) — /3(a) parameterisation can also be used to 
analyse the gravitational tests in the solar system. For the 
chameleon models, i.e. the first example, and evaluating 
Eq. (11) in the galactic vacuum, we find that 



topi 



9/3 Q n r 



H 2 



da- 



9(a) 



where aa ~ 10 -2 is the scale factor when the matter den- 
sity in the cosmological background equals the galactic 
density pc ~ 10 6 p c . Here 4>a is the value of the scalar 
field in the galaxy and <f> c the field inside ordinary mat- 
ter. Local tests are satisfied when the thin shell effect is 
at play where we define the thin shell factor 
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4>g - 4>c 

6to p1 ^ c $ g 
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and the modification of gravity in the solar system felt 
by a satellite such as the Cassini probe has a strength 
2/?g/?c^tj^ in which (3 G is the coupling in the galactic 
vacuum, $ Q is the value of the solar Newtonian poten- 
tial ($© ~ 1CT 6 ) and f3 c is the coupling inside a dense 
body The overall result should be less than 10 -5 to com- 
ply with the Cassini bound in the solar system [13]. This 
condition is independent of [3 C and reads 



Lucul CiuiMmiiil* 



G 



da 



9(a) 
a 4 / 2 (a) 



< 10" 



9a 



<l>: 



The integral / = J a "° da a ff2( a \ is potentially divergent 
for small values of a; n j ~ 10~ 10 . Hence we must impose 
that f(a) 2 /g(a) compensates the 1/a 4 divergence in the 
integrand. Typically, we can parameterise 



/(«) 



9(a) 



We must then impose 2r — s > 3, which leads to I ss 
The case s = and r > 3/2 has already be seen 



2r-s-3 1 

to correspond to inverse power law chameleons and large 
curvature f(R) models. Adopting j3c w 10 2s /3o, we find 
that there is an interplay between /3q and mo/Ho: 



< io 



4r-4s-12 



(2r- 



3) 



In the above we have assumed that galaxies have a thin 
shell to minimise the disruption of their dynamics, al- 
though the necessity of this condition should be ascer- 
tained using N-body simulations [14]. Enforcing the thin- 
shell condition imposes 



\4>g-M < 6/3 $G 

where the galactic Newtonian potential is $g 
and 



io- 



— — 9/?o^mO — % 

mpi mi 



da- 



9(a) 
« 4 / 2 («) 



implying that 



^>10 3 



(23) 



where </>o is the cosmological value of <fr now. This con- 
dition is independent of /3q and means that any screened 
modified gravity model will have effects on Mpc scales 
only. 

Strong constraints can also be obtained from labora- 
tory experiments. Using the fact that the initial matter 
density at Zj n j ~ 10 10 corresponds to the matter density 
inside laboratory test bodies, we have mi a b ~ 10 10r mo. 
Gravity is not modified provided tests bodies have a thin 
shell, 



\4> 



lab 



< 6/3labm P1 $ 



10 s - 



io"- 




lab 



FIG. 1: The constraints on mo /Ho as a function of r for /3o = 
l/v6 and s = 0. Valid models must be above the red (solar 
system), mauve (cavity) green (m > H), light blue (md > 1), 
light red (fi) and brown (galaxy) lines. The blue line gives the 
detectability of effects on the CMB by the Planck satellite. 
The strongest constraints are the cavity and galactic bounds 
for small and large r respectively. Models with r > 3 satisfy 
the constraints and can lead to a modified gravity regime on 
large scales. 



where $i a b ~ 10~ 27 for typical test bodies in cavity ex- 
periments of size R, 0i a b = 4>(auh) where m(ai a b) ~ 1/R 
. A weaker condition is that mi a bd 3> 1 where d is the 
size of the test body implying that mod ^> 10 _10r . Both 
constraints can be found in figure 1. Finally, the scalar 
field mass is larger than the Hubble rate since a = a- m i 
provided r>2 — ln^|^ / lna ini . 

For the dilaton models, if the coupling now is of or- 
der unity, according to Eq. (22), (3(4> c ) ^ IO -2 5 can 
be achieved provided that A-i > and that the time 
variation of m(a) is slow and does not compensate the 
1/a 4 divergence in the integrand. In this situation, the 
coupling function (3 converges exponentially fast towards 
zero: this is the Damour-Polyakov mechanism [15]. Al- 
ternatively, a smooth variation of the coupling to mat- 
ter and therefore interesting consequences for large scale 
structures are achieved when the mass of the scalar field 
compensates exactly the 1/a 4 . This is obtained for mod- 
els with to 2 (a) = 3A2-ff 2 (a) where A2 3> 1 here. Indeed, 
-ff(a) ~ 1/a 2 in the radiation era implying that the time 
variation of (3 until the matter-radiation equality is very 
small. In the matter era H (a) ~ 1/a 3 / 2 implying a power 
law variation of /3 with a. These models have been ex- 
tensively studied and correspond to the environmentally 
dependent dilatons [10] where effects of modified gravity 
are also only effective at the Mpc scale. 

The local constraints give that mo /Ho > 10 3 which 
implies that astrophysical effects of modified gravity can 
only occur on Mpc scales. This fact can be easily seen by 
stydying the growth of structures in the linear regime. In 



5 



the matter dominated era, the density contrast of cold 
dark matter (CDM) evolves according to 



5"+U8' c -lu 2 P — ( 1 + - 

2 Pc + P-, + Pb\ 1 



2/3 2 



0. (24) 



Hence inside the Compton radius k > am (a), perturba- 
tions grow anomalously. Modified gravity can be effective 
if in the recent past of the Universe the Compton radius 
is of the order of 10 Mpc, i.e. mo/H > 10 3 . We have 
seen that this is a consequence of local constraints for 
screened modified gravity models. On the other hand, on 
such scales non-linear effects cannot be neglected. The 
strength of the m(a) — /3(a) paramcterisation is that the 
non-linear regime can be easily described. Matter cluster- 
ing on galactic and cluster scales is an important probe of 
modified gravity. On Mpc scales, the nonlinearity in both 
the structure formation process and the dynamics of the 
scalar field requires full numerical simulations [18, 19]. 
The /3(a), m(a) parameterisation can completely specify 
the nonlinear dynamics of <fi. To see this, note that in the 
quasi-static limit the scalar field is governed by 

V 2 = -{a^p m ~ a^pm) + V^(0) - V,,f,(4>), (25) 

where the overbar means the background value. With 
4>{a) reconstructed from /3(a) and m(a), and V^ia) from 
/3(a) and p(a), one can easily obtain V^(0) analytically 
or numerically, and this can be used to solve the quasi- 
static dynamics numerically. An advantage is that tempo- 
ral functions m(a), /3(a) completely specify the dynamics 
of (f>, in particular its spatial configuration, and there is 
no need for a fc-space parametrisation. 

On the other hand, linear scales could also be of in- 
terest for deciphering modifying gravity if, for scales 
entering the horizon before radiation-matter equality, 
the anomalous growth of perturbations plays a role on 
the peak structure of the cosmic microwave background 
(CMB). This happens provided the ratio 



(i + '^y/ 



is of order one [16] . The scalar field mass at last scattering 
is given by acMB m CMB ~ 10 3r_3 mo, hence we find that 



/3cmb 



fccMB^CMB ^ -|^3-3r+3s ^CMB PqHq 



(!CMB"*CMB H m Q 

where fccMB characterises the scale of the horizon at the 
last scattering. The modified gravity effects can be seen 
on the CMB provided 10 2 < u « l/(2/3 2 MB ) < 3 ■ 10 3 
[17], giving constraints on the model parameters through 
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2uj kcMB 



— io 3( - 1+r ~ s) 



in which typically we have Hg/kcMB ~ O(0. 1). The pre- 
vious condition on the observability of modified grav- 
ity by the Planck satellite is compatible with the solar 



system constraint provided r — s < —7 and therefore 
PoHo/ttiq < 10~ 23 . This implies that for reasonable val- 
ues of /3o, mo/Ho would be so large that no effect of 
modified gravity on large scale structures (LSS) would 
be present. Hence, for models with a power law depen- 
dence of both the mass and the coupling to matter, effects 
on both the CMB and LSS are not compatible. 

The scalar field also has an effect on gauge couplings 
and particle masses. The fermion masses are given by 
m^(</>) = A(<fi)m F where m F is the bare mass in the La- 
grangian. Meanwhile, quantum effects such as the pres- 
ence of heavy fermions lead to the coupling of to pho- 
tons 

Sgaugc = J d 4 X^B F (<t>)F^F^, 

where g is the bare coupling constant and 



m Pi 



(26) 



Depending on the model, the coefficients /3 and /3 7 can be 
related. Here we will consider them to be free parameters. 

The scalar coupling to the electromagnetic field could 
lead to a dependence of the fine structure constant on <f> 
as 



implying that 



1 

Q 



Q 
Q 



1 

ao 



Bf {<!>), 



-/3 7 K4</>. 



Using the evolution equation we find that 



« -9/3 7 /3fi m - 
Ha m 



(27) 



(28) 



(29) 



Hence the negative variation of the fine structure con- 
stant in one Hubble time is related to the small ratio 
S < 1 and the couplings of (j> to matter and photons. 
The best experimental bound on the variation of a now 
comes from Aluminium and Mercury single-ion clocks: 
(-1.6±2.3)-10- 17 



™ o 



yr 1 . Taking H^ 1 - 1.5-10 10 yr, 



we get the conservative bound 



Ha 



< 2 • 10" 7 . 



As a result, the experimental bounds on the time vari- 
ation of a lead to constraints on /3o/3 7 o as /3o/3 7 o < 

0.8 • 10" 7 ^|. For models with /3 = 0(1), fi m - 0.25 

and mo /Hq w 10 3 where effects on LSS are present, this 
is a tighter bound than present experimental ones [24] 



< 0.1 



(30) 
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Fundamental fcrmions such as the electrons have a uni- 
versal mass dependence mp = A((f>)m F , implying that 



mp 
Hnip 



H 2 

9(3 2 n m —. 



HI 



(31) 



Similarly, nucleons such as the proton have a mass given 
by the phenomenological formula 

m v = CqcdAqcd + b u m u + b d m d + C p a, 

where Aqcd is the QCD scale, b u + bd ~ 6, b u — bd^ 0.5, 
Cqcd ~ 5.2, m° - 5MeV, m° d ~ lOMeV and C p a ~ 
0.62MeV. Of course the main source of uncertainty here 
follows from the lack of our knowledge about the coupling 
of scalars to gluon. This leads to weaker bounds than the 
ones coming from the time variations of a. Assuming 
conservatively that Aqcd is scalar independent 1 , we get 



Hm, 



H 2 n u m° u 



bdm d 



(32) 



It is particularly important to study the variation of \x = 
— from which we find that its time variation is positive 
for modified gravity models: 



it m 



m p 



(33) 



The current experimental constraint is = (—3.8 ± 
5.6)10 _14 yr _1 which yields the upper bound on j3g: 



(34) 



For /3 = 0(1), this entails that m /H > 10 25 . In Figure 
1, we summarise all the constraints when s = and /3o = 
l/y/E (corresponding to inverse power law chameleons 
and large curvature f(R) models) and show that models 
with r > 3 (large curvature f(R) models) are compatible 
with effects of modified gravity on LSS. 

Modified gravity models compatible with local experi- 
ments must inherit some mechanism to suppress a poten- 
tial fifth force. This makes these models very nonlinear. 
We have shown that a large class of such models can be 
fully parameterised by only two temporal functions: the 
scalar field mass and the coupling to matter allowing one 
to reconstruct the full action using (8) and (10). An im- 
portant implication is that instead of studying individual 



1 If quantum effects lead to a coupling of <j> to gluons with 
B(<f>) ~ 1 + (3 g K4(j> like for photons, then Aqcd oc (1 + 



/3 S K40 [25], where a s (M z ) 



0.12 



is the running strong coupling at the energy scale of the weak Z 
boson mass. The time variation of Aqcd can be much faster than 
that of a, leading to stronger constraints. For a recent analysis, 
see [26] 



models, one can focus on the effects entailed by the choice 
of these two functions. We have shown that many known 
models can be reconstructed by using a simple power-law 
form for these two functions. Solar system tests, labora- 
tory experiments and the variation of fundamental cou- 
plings and masses can be easily studied and lead to strong 
constraints on models. In cosmology, the CMB is rather 
weakly influenced by the fifth force, while on Mpc scales 
the LSS can be largely affected. On such scales the non- 
linearity of both the structure formation and the models 
is important, and N-hody simulations are needed to fully 
understand the model behaviour. This work is ongoing. 

We would like to thank Justin Khoury and Jean- 
Philippe Uzan for useful discussions. 
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